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Abstract 

In this article, we derive the stochastic master equations corresponding to the sta- 
tistical model of a heat bath. These stochastic differential equations are obtained as 
continuous time limits of discrete models of quantum repeated measurements. Physically, 
they describe the evolution of a small system in contact with a heat bath undergoing 
continuous measurement. The equations obtained in the present work are qualitatively 
different from the ones derived in [B], where the Gibbs model of heat bath has been studied. 
It is shown that the statistical model of a heat bath provides clear physical interpretation 
in terms of emissions and absorptions of photons. Our approach yields models of random 
environment and unravelings of stochastic master equations. The equations are rigorously 
obtained as solutions of martingale problems using the convergence of Markov generators. 

Introduction 

The theory of Quantum Trajectories consists in studying the evolution of the state of an 
open quantum system undergoing continuous indirect measurement. The most basic physical 
setting consists of a small system, which is the open system, in contact with an environment. 
Usually, in quantum optics and quantum communication, the measurement is indirectly per- 
formed on the environment [BJ [13l [TU1 [28l l39l I40j . In this framework, the reduced time 
evolution of the small system, obtained by tracing over the degrees of freedom of the environ- 
ment, is described by stochastic differential equations called stochastic Schrddinger equations 
or stochastic Master equations. The solutions of these equations are called Continuous Quan- 
tum Trajectories. In the literature, two generic types of equations are usually considered 

1. Diffusive equations 

d Pt = C( Pt )dt + [Cp t + Pt C* - Tr [p t (C + C*)] p t )dW t , (1) 
where (Wt)t>o is a one dimensional Brownian motion. 
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2. Jump equations 

dp t = C(p t )dt + (^^—^ - p^j (dN t - Tr [C pt C*) dt) , (2) 

where (Nt)t>o is a counting process with stochastic intensity t \—* f Tr [C p s C*~\ds. 

Physically, equation (TjQ) describes photon detection models called heterodyne or homodyne 
detection [7J EJ G22 HO] . The equation ((2J) relates direct photon detection model [391 HO] • 
The driving noise depends then on the type of measurement. Mathematically, a rigorous 
approach for justifying these equations is based on the theory of Quantum Stochastic Calculus 
[SI ESI El IM] • In such a physical setup, the action of the environment (described usually by a 
Fock space) on the small system is modeled by quantum noises [HO [21]. The evolution is then 
described by the so-called Quantum Stochastic Differential Equations [1, 2, 24, 20J. Next, by 
using the quantum filtering \16\ \18\ IT9] technique, one can derive the stochastic Schrodinger 
equations by taking into account the indirect observations. Another approach, not directly 
connected with quantum stochastic calculus, consists in using instrumental operator process 
and notion of a posteriori state [71 151 ITUl ITU fT2l [22] . 

In this work, we shall use a different approach, introduced recently by the second author 
in |25 | I26 | 127]. This discrete-time model of indirect measurement, called Quantum Repeated 
Measurements is based on the model of Quantum Repeated Interactions [3 HI [S] introduced 
by S. Attal and Y. Pautrat. The setup is the following: a small system TC is in contact with 
an infinite chain, (^j^^fc, of identical and independent quantum systems, that is £). = £ 
for all k. The elements of the chain interact with the small system, one after the other, 
each interaction having a duration r > 0. After each interaction, a quantum measurement is 
performed on the element of the chain that has just been in contact with the small system. 
Each measurement involves a random perturbation of the state of the small system, the 
randomness being given by the outcome of the corresponding quantum measurement. The 
complete evolution of the state of the small system is described by a Markov chain depending 
on the time parameter r. This Markov chain is called a Discrete Quantum Trajectory. By 
rescaling the intensity of the interaction between the small system and the elements of the 
chain in terms of r, it has been shown in |25| [2^] that the solutions of equations 111 121) can be 
obtained as limits of the discrete quantum trajectories when the time step r goes to zero. 

In |25| HH] j the author investigated the case when the reference state of each element of the 
chain is the ground state (this corresponds also to models at zero temperature). This setup 
was generalized in [6], where Gibbs states with positive temperature were considered and the 
corresponding equations were derived. In the present work, we go beyond this generalization 
and study the statistical model for the temperature state of the chain. More precisely, the 
initial state of the elements of the chain is a statistical mixture of ground and excited states. 
It is important to notice that both the Gibbs model as well as the ground state model are 
deterministic. Let us stress that, in the case where no measurement is performed after each 
interaction, both the Gibbs and the statistical model give rise to the same deterministic limit 
evolution. This limit behavior confirms the idea that a mixed quantum state and a proba- 
bilistic mixture of pure states represent the same physical reality. Quite surprisingly, we show 
that, when adding measurement, the limit stochastic differential equations are of different 
nature: for the Gibbs model the only possible limit evolutions are deterministic or diffusive, 
whereas for the statistical model jump evolutions becomes a possibility. Furthermore, the 
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Gibbs model limit equations involve at most one random noise, whereas two driving noises 
may appear at the limit when considering in the statistical model. 

The article is structured as follows. In Section HJ we introduce the different discrete 
models of quantum repeated interactions and measurements. In our approach, we present 
the statistical model of the thermal state as the result of a quantum measurement applied to 
each element of the chain before each interaction. Next, we describe the random evolution of 
the open system by deriving discrete stochastic equations. In Section [2j we investigate the 
continuous time models obtained as limits of the discrete models when the time-step parameter 
goes to zero. We remind the results of [6] related to the thermal Gibbs model and we describe 
the new continuous models related to the thermal statistical model. Section [3] is devoted 
to the analysis of the different models. The qualitative differences between the continuous 
time evolutions are illustrated by concrete examples. Within these examples, it is shown that 
the statistical approach provides clear physical interpretations which cannot be reach when 
considering the Gibbs model. We show that model at zero temperature (each element of the 
chain is at the ground state) can be recovered by the statistical model; however, this is not 
possible with the Gibbs model. Moreover, we show that considering the statistical model 
allows to obtain unravelings of heat master equations with a measurement interpretation. 
Section 0] contains the proofs of the convergence of the discrete time model to the continuous 
model. Such results are based on Markov chain approximation techniques using the notion 
of convergence of Markov generators and martingale problems. 

1 Quantum Repeated Interactions and Discrete Quantum Tra- 
jectories 

In this section we present the mathematical model of quantum repeated measurements. In 
the first subsection we briefly recall the model of quantum repeated interactions [I] and in the 
second subsection we describe three different situations of indirect quantum measurements, 
in which environment particles are measured before and/or after each interaction. Discrete 
evolution equations are obtained in each case. 

1.1 Quantum Repeated Interactions Model without Measurement 

Let us introduce here the mathematical framework of quantum repeated interactions. We 
consider a small system TC in contact with an infinite chain of identical and independent 
quantum systems. Each piece of the chain is represented by a Hilbert space £. Each copy of 
£ interacts, one after the other, with the small system TC during a time r. Note that all the 
Hilbert spaces we consider are complex and finite dimensional. 

We start with the simpler task of describing a single interaction between the small system 
TC and one piece of the environment £. Let p denote the state of TC and let a be the state of 
£. States are a positive self-adjoint operators of trace one; in Quantum Information Theory 
they are also called density matrices. The coupled system is described by the tensor product 
TC® £ and the initial state is in a product form p ® a. The evolution of the coupled system 
is given by a total Hamiltonian acting on TC <S> £ 

H tot = H ®I + I®H + flint, 

where the operators flo and fl are the free Hamiltonians of the systems TC and £ respectively, 
and the operator H mt is the interaction Hamiltonian. The operator fl to t gives rise to an 
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unitary operator 



U = exp(-«7\fftot), 



where r represents the time of interaction. After the interaction, in the Schrodinger picture, 
the final state of the coupled system is 

p = U{p®a)U*. 

In order to describe all the repeated interactions, we need to describe an infinite number 
of quantum systems. The Hilbert space of all possible states is given by the countable tensor 
product 

oo 

r = H®§§e k = H®®, 

k=l 

where £ & ~ £ for all k > 1. If {eo, e±, . . . , ex} denotes an orthonormal basis of £ ~ C K+1 , 
the orthonormal basis of $ = (^feLi £j, is constructed with respect to the stabilizing sequence 
e o N ( we snan n °t develop the explicit construction of the countable tensor product since we 
do not need it in the rest of the paper; we refer the interested reader to [1] for the complete 
details) . 

Let us now describe the interaction between H and the k-th piece of environment £/., 
from the point of view of the global Hilbert space T. The quantum interaction is given by 
an unitary operator Uk which acts like the operator U on the tensor product 7i®£k and like 
the identity operator on the rest of the space T. In the Schrodinger picture, a state rj of T 
evolves as a closed system, by unitary conjugation 

rj^U k r,Ul 

Therefore, the whole procedure up to time k can be described by an unitary operator V k 
defined recursively by 

/ Vfc+i = Uk+iVk 

\Vo = I {6) 

In more concrete terms we consider the initial state /J, = p (g> <r for the small system 

coupled with the chain (notice that all the elements of the chain are initially in the same state 
= cr )- After k interactions, the reference state is given by 

fJ-k = VuiiVk- 

Since we are interested only in the evolution of the small system TC, we discard the 
environment 3>. The reduced dynamics of the small system is then given by the partial trace 
on the degrees of freedom of the environment. If a denotes a state on T, we denote by Tr^[a] 
the partial trace of a on H with respect to the environment space $ = ^^=1 ^k- We recall 
the definition of the partial trace operation. 

Definition-Theorem 1 Let TC and tC be two Hilbert spaces. For all state a onH^K,, there 
exists a unique state on TC denoted by Tr/c [a] which satisfies 

Tr [Tr/c[a]X] = Tr[a(X ® I*;)], 

for all X £ B(7i). The state Tr^;[a] is called the partial trace of a onTC with respect to 1C. 
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With this notation, the evolution of the state of the small system is given by 

p k = Tr* [nk\- ( 4 ) 
The reduced dynamics of (p k ) is entirely described by the following proposition [U [25] . 

Proposition 1 The sequence of states (pk)k defined in equation satisfies the recurrence 
relation 

Pk+1 = Tr £ [U(p k ®<r)U*]. 

Furthermore, the application 

L : B(H) -» B(H) 

X h+ Tr e [U(X ® a)U*\ 

defines a trace preserving completely positive map (or a quantum channel^ and the state of 
the small system after k interactions is given by 

Pk = L k (p ). (5) 
1.2 Quantum Repeated Interactions with Measurement 

In this section we introduce Quantum Measurement in the model of quantum repeated in- 
teractions and we show how equation ([5]) is modified by the different observations. We shall 
study three different situations of indirect measurement, as follows: 

1. The first model concerns "quantum repeated measurements" before each interaction. 
It means that we perform a measurement of an observable on each copy of £ before the 
interaction with H. We call such a setup "Random Environment" (we shall explain the 
terminology choice later on). 

2. The second model concerns "quantum repeated measurements" after each interaction. 
It means that we perform a measurement of an observable on each copy of £ after the 
interaction with H. We call such a setup "Usual Indirect Quantum Measurement". 

3. The third setup is a combination of the two previous models. Two quantum measure- 
ments (of possibly different observables) are performed on each copy of £, one before 
and one after each interaction with H. Such a setup is called "Indirect Quantum 
Measurement in Random Environment" 

In all the cases, the measurement is called indirect because the small system is not directly 
observed, the measurement being performed on an auxiliary system (an element of the chain) 
which interacted previously with the system. The main purpose of this work is to study and 
analyze the three different limit behaviors obtained when the interaction time r goes to zero 
(see Section 2) . Let us mention that the second setup has been studied in detail in [251 EH EZ] • 
We chose to describe in great detail the more general case of the third model, since the other 
two models can be easily recovered from the third one, by choosing to measure the trivial 
observable I. 
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1.2.1 Indirect Quantum Measurement in Random Environment 

In order to make the computations more easy to follow, we shall focus on the case where the 
environment is a chain of qubits (two-dimensional quantum systems). Mathematically, this 
is to say that £ = C 2 . 

Let us start by making more precise the physical model for one copy of £. To this end, 
we consider {eo,ei} an orthonormal basis of £, which diagonalizes the Hamiltonian 



H = 



7o 
71 



where we suppose that 70 < 71. The reference state a of the environment corresponds to a 
Gibbs thermal state at positive temperature, that is 

e -pH 1 

a = T^FPry with fi = kt> (6) 

where T corresponds to a finite strictly positive temperature and if is a constant. In the 
basis {eo,ei}, a is diagonal 

o =p|e )(e | + (1 -p)|ei)(ei|, 

with 

e -/370 

P ■ 



e -Plo + e -Pn " 

Notice that since [3 > 0, we have < p < 1. 

We are now in position to describe the measurement before the interaction. We consider 
a diagonal observable A of £ of the form 

A = A |e )(e | + Ai|ei)(ei|. 

The extension of the observable A to an observable of H®£ is 1(8) A. According to the axioms 
of Quantum Mechanics, the outcome of the measurement of the observable I<S)A is an element 
of its spectrum, the result being random. If the initial state (before the interaction) is p<g> a, 
we shall observe the eigenvalue Aj with probability 

P[A; is observed] = Tr [(p ® a) I (8 Pi] = Tv[aPi\, i = 0, 1 

where Pi = |ej)(ej| are the eigenprojectors of A. It is straightforward to see that in this case 

P[Ao is observed] = p = 1 — P[Ai is observed]. 

Furthermore, according to the wave packet reduction principle, if the eigenvalue Aj is observed, 
the initial state p (8 o is modified and becomes 

^ Tr[( P 0a) P ^Tv[aP^ UJ 

This defines naturally a random variable p 1 valued in the set of states on ?i®£. More precisely, 
the state p l takes the value p\ = p (8 |eo)(eo| with probability Tr [(p (8 <r) |eo)(eo|] = p and 
the value p\ = p (8 |ei)(ei| with probability 1 — p. 
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Remark 1 Since both the initial state of the system and the observable measured have product 
form, only the state of £ is modified by the measurement before the interaction. Instead of 
describing the evolution of the coupled system, we could have considered that the state of £ is a 
random variable a\ where o~\ is either |eo)(eo| with probability p either |ei)(ei| with probability 
1 — p. This is the statistical model for a thermal state and its random character justifies the 
name "Random environment". In conclusion, we could have replaced from the start the setup 
(Gibbs state + Quantum measurement) with the probabilistic setup Random environment, 
the results being identical. We shall give more details and comments on this point of view in 
the Section^ 

We now move on to describe the second measurement, which is performed after the inter- 
action. In this case we consider an arbitrary (not necessarily diagonal in the basis {eo,ei}) 
observable B of £ which admits a spectral decomposition 

B = a Qo + aiQi, 

where Qj corresponds to the eigenprojector associated with the eigenvalue ay. Let /x 1 be the 
random state after the first measurement. After the interaction, the state on Ti ® £ is 

r,\ = U^U\ £ = 0,1. 

Now, assuming that the measurement of the observable A (before the interaction) has given 
the result Aj, the probability of observing the eigenvalue ay of B is given by 

P[ay is observed] = TV [nj I <g)Qj] . 

and the state after the measurement becomes 

0i _ I ®Qj yj I ®Qj 
i>j ~ Tr [ril I ®Q-] ' 

The random state 6 l (which takes one of the values 6jj) on7i(g>£ describes the random result 
of the two indirect measurements which were performed before and after the interaction. 

Having described the interaction between the small system TC and one copy of £ , we look 
now at the repeated procedure on the whole system T. The probability space underlying the 
outcomes of the repeated quantum measurements before and after each interaction is given 
by f2 = (£a x S#) n *, where = {0,1} corresponds to the index of the eigenvalues of 
the observable A and = {0, 1} for the ones of B. On f2, we consider the usual cylinder 
o"-algebra A generated by the cylinder sets 

A (h,ji),...,(i k ,jk) = e x £ B ) N \ui = ii,...,Uk = ik,<Pi =ji, ■ ■ ■ ,<Pk = 3k}- 

Now, we shall define a probability measure describing the results of the repeated quantum 
measurements. To this end, we introduce the following notation. For an operator Z on £j, 
we note Z^> the extension of Z as an operator on T, which acts as Z on the j-th. copy of £ 
and as the identity on 7i and on the other copies of £: 

3-1 

Z® = I®(g)I®Z® (g) I. 

p=i p>j+i 
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Furthermore, for all k > 1 and {(h,ji), • • • , (ijfc, jk)} S (£4 x £b) , we put 





where Pi and are the respective eigenprojectors of A and B and = p® ^^=1 with 
<Tfc = a = p|eo)(eo| + (1 — p)|ei)(ei| for all k G N*, is the initial state on F. Notice that the 
products in the previous equation need not to be ordered, since two operators and Y^' 
commute whenever i ^ j. In the same vein, the following important commutation relation 

QMu k pW... Q M UlP U 

shows that the operator p k ({h,ji), ■ ■ ■ , (ik,jk)) in Eq. (jHJ) is actually the non normalized state 
of the global system after the observation of eigenvalues Aj x , . . . , Aj fc for k first measurements 
of A and a^, . . . , ctj k for the k first measurements of the observable B. 

We have now all the elements needed to define a probability measure on the cylinder 
algebra A by 

P [ A (ii,ii),...,(i fc ,i fe )] = Tr[/i fe ((n, ji), . . . , (i fc , j fc ))]. 

This probability measure satisfies the Kolmogorov Consistency Criterion, hence we can extend 
it to the whole u-algebra A to the unique probability measure P with these finite dimensional 
marginals. 

The global random evolution on V is then described by the random sequence (p k ) 
p k : Q — B(T) 

I x ~ I \ A((^l, tpi), • • • , (Vk,<Pk)) 



Tr[/t fc ((o;i,^i), . . . , (wfc,^))] 

This random sequence describes the random modification involved by the result of measure- 
ment before and after the interactions. In order to recover the measurement setup only before 
or only after the interactions, one has just to delete the projector P} or Qf, in equation ([8]). 
The reduced evolution of the small system is obtained by the partial trace operation: 

Pk (uj, Lp) = Tr$ [p k (u, (p)] (9) 

for all (cj, ip) S Q and all fc£N*. The random sequence (p k )k>i is called a Discrete Quantum 
Trajectory. It describes the random modification of the small system undergoing the sequence 
of successive measurements. 

Remark 2 The dynamics of the sequence of states pk can be seen as a random walk in 
random environment dynamics in the following way. Assume that all the elements of the 
chain are measured before the first interaction; the results of this procedure define a random 
environment in which the small system will evolve. All the randomness coming from the 
measurement before each interaction is now contained in the environment uj. Given a fixed 
value of the environment lo , the small system interacts repeatedly with the chain (whose states 
depend on lo) and the random results of the repeated measurement of the second observable 
B are encoded in ip. In this way, the global evolution of p can be seen as a random walk 
(where random modifications of the states are due to the second measurement) in a random 
environment (generated by the measurements before each interaction). 
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1.2.2 Discrete Evolution Equations 

In this section, using the Markov property of the discrete quantum trajectories, we obtain 
discrete evolution equations which are random perturbation of the Master equation ([5]) given 
in Proposition [H The Markov property of the random sequence (p k )k is expressed as follows. 

Proposition 2 The random sequence of states (pk)k on T~t defined by the formula (0 is a 
Markov chain on (fl, A, P). More precisely, we have the following random evolution equation 

where 

Gijip) = Tr £ [(I®Qj) U {l®Pi {p®a) I®Pi) U* {I®Qj)} 
and 1^(0;,^) = lij(u k+1 ,ip k+1 ) = li(u k+1 )lj(ip k+1 ) for all (u,ip) G (£4 x £ S ) N *. 



The equation (I10p is called a Discrete Stochastic Master Equation. In order to make more 
explicit the equation (fTO]) and to compute the partial trace, we introduce a suitable basis for 
H (g) £, which is {eo eo, ei (8) eo, eo (8> ei, ei (8) ei}. In this basis, the unitary operator £/ can 
be written in block format in the following way 



U 



Loo Loi 
Lio Lu 



where L« are operators in A^2(C). We shall treat two different situations, depending on the 
form of the observable B that is being measured after each interaction. On one hand we 
consider the case where the observable B of E is diagonal in the basis (eo,ei) and on the 
other hand we consider the case where B is non diagonal. 

Let us start with the case where the observable B is diagonal in the basis {eo, ei}, that is 
B = ao|eo)(eo| + «i|ei)(ei|. In this case, equation (fTQ|) becomes 

p k+1 (u,<p) = 

L op k {uJ,ip)L^ L 1Q p k {u,y)L\ Q 



Tr[L o/5fc(^, <^)-Loo] ' Tr[Li pfc(w,(/?)-Li ] 

Lpipfc(a^y)L^ Lu p k {u>, if) L^ 

Tr[L i p fc (w, v?) LjJJ Tr[Ln p fe (u;, y?) L^] 

(11) 

Usually, a stochastic Master equation appears as a random perturbation of the Master 
equation (see equations ([HE]) in the Introduction). Moreover, the noises driving the equations 
are centered, that is of zero mean (this is the case of the Brownian motion and the counting 
process compensated with the stochastic intensity in equations ((HE]))- In order to obtain a 
similar description in the discrete case, we introduce the following random variables 

X k {u>,(p) = , ken. (12) 

y/p{l-p) 
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Now, we rewrite equation (jlip in terms of the random variables X k , loi lio : 
Pk+i = (p{LooPkLoo + L 10 p k L* ) + (1 -p)(L m p k L* m + L u p k L* n ) )l 



- vV(l -p) (Loo Pk Lq + L 10 p k L* 10 ) + y/p(l -p) (L u p k L* n + L 01 p k L^) )X k+1 

I L ooPk L* L wPk L* \ 

+ { ~ TY[WZ4l + nL 10 p k L to ] ) (l01 " P ^ W*o]) 

+ ( ~ Jl' 1PkL % ] + ^T^%t ) (lio " (1 " P) Tr[L 0lPfc L 01 ]). (13) 

It is important to stress out that the last three terms in the previous equation have mean 
zero: 

E[X k ] =E[l i -pTr[L 10 p fc L 10 ]] = E[l 10 - (1 - p) Tr[L ip k L 01 ]] = 0. 

Moreover, recall that the discrete evolution of Proposition Q3 without measurement, is given 
by 

Pk+i = L(p k ) = (p(L QO p k Lq + L w p k L* ) + (1 -p)(L 01 p k V m + L u p k L\^j. (14) 



As a consequence, the discrete stochastic master equation (jl3|) is written as a perturbation 
of the discrete Master equation (fill . 

Remark 3 In this expression, one can see that the random variable X k depends only on the 
outcome of the measurement before the interaction (we sum over the two possible results of the 
measurement after the interaction). In other words, it means that the random variables X k 
describe essentially the perturbation of the measurement before the interaction. On the other 
hand, the random variables loi and conditionally on the result of the first measurement, 
describe the perturbation involved by the measurement after the interaction. Hence, each term 
of the equation t!3\) that is linked with either X k , loi or lio expresses how the deterministic 
part is modified by the results of the different measurements. 

We now analyze the second case, where the observable B is non-diagonal in the basis 
{eo, ei}. We write B = aoQo + «iQl> where the eigenprojectors Qi are written in the {eo, e{\ 
basis Qi = (o%i)o<k,l<l- I n this case, the operators appearing in equation (fTU|) are given by 

Goi(p) = Qoo l oopLq + ql L 00 pLl + q^LiopL^Q + q\±Li pL* 
Oli(p) = Qoo L oipLoi + qloLoxpL*! + q^LnpL^ + q'^LnpL^. 

As before, in order to obtain the expression of the discrete Master equation as a perturbation 
of the deterministic Master equation, we introduce the following random variables 



x. 



k+1 



Y° 
r k+l 



l^-pTVfSoiM 



pTr[g Q1 (p k )](l-pTr[g Q1 (p k )]) 
lt +1 -(l-p)Tr[g 10 (p fc ) 

(1 - p) tv [g w (p k )] (1 - (1 - p) tv [g w (Pk)] ) 



n 1 1 = i?o +1 -(i-p)Tv[g 10 (p fc )] (i5) 
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In terms of these random variables, we get 



fc+i 



Pk+i = L(p k )l 

+ ( - VpO- -P){Goo{Pk) + Gm{pk)) + VpO- -p)[Gu(pk) + Gio(Pk))) Xk+i 

« -'™»<*»> ( - + Y ° 

Wd - rt WW -d-,) »)]) ( - + ^gL) nv, 

(16) 

Remark 4 ^4s it was the case in equation $13\) . the discrete random variables X k and Y k , 
i = 0,1 are centered. As before, the variables X k represent the perturbation produced by the 
measurement before the interaction and, given the result of this measurement, the variables 
Y k % describe the perturbation generated by the measurement of the second observable. The 
particular choices made for Y k will be justified when we shall consider the continuous models. 
They will appear as discrete analogs of the noises which drive the continuous stochastic Master 
equations (Wt and N t in equations 

The above general framework concerns the combination of the two measurements, one 
before and one after each interaction. Let us present the corresponding equations when only 
one type of measurement (before or after each interaction) is performed. 

We start by looking at the case where a measurement is only performed before the inter- 
action (we called this kind of setup "Random environment"). Since measuring an observable 
on an element of the chain (which has not yet interacted) does not alter the state of the 
little system 7i, only the reference state of each copy of £ is random. The completely positive 
evolution operators describing the two possibilities for the state after the interaction are given 
by 

Ki{p) = Ti £ [U(p®\e l ){e i \)W) (17) 
= Gio(p) + Ga(p), (18) 

for i = 0, 1. Let 1^ be the random variable which is equal to 1 if we observe the eigenvalue 
Aj at the k-th step, and otherwise. We can describe the evolution of the small system H by 
the following equation 

Pk+1 = n (p k )i k Q +1 +n 1 ( Pk )i k 1 +1 . 

As before, we introduce 

vp(i -p) 

With this notation, the evolution equation becomes 



Pfc+i = L(p k ) + ( - y/p(l -p) (Lqo Pk Lq + L 10 p k L\ ) 



+VW - P)(Ln Pk L* n + Loi Pk L* m ))X k+1 . (19) 

The opposite case, where a measurement is only performed after the interaction, is treated 
in great detail in [251 ES] when p = (ground states) and in [6] for < p < 1 (Gibbs states). 
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Let us recall briefly the main steps needed to obtain the appropriate equations. Consider the 
observable B = a.\Q\ + 012Q21 with Qi = (?jy)o<fc,J<i- The two possible non normalized states 
on TC that can be obtained after the measurement are defined via the action of the operators 

Fi{p) = TT £ [I®Q i U(jp\e Q ){e Q \ + {l-p)\ex){e x \)U*I®Q i } (20) 
= pGoi(p) + {l-p)Gu(p), (21) 
for i = 0,1. The discrete evolution equation is then given by 

Again, we introduce the random variables defined by 

lf^-Trf^)] 



Xk+i 



VTr^P^Tr^O^)] 
In terms of these centered random variables, we get 




Pk+1 = H Pk )l +\-J ™ [^| MPk) + ^ [g^j MPk) j X h+1 . (23) 
2 Continuous Time Models of Quantum Trajectories 

In this section, we present the continuous versions of the discrete equations (fT3l [TBI [191 123|) . 
We start by introducing asymptotic assumptions for the interaction unitaries in terms of the 
time parameter t. Next, we implement these assumptions in the different equations <\13\ [To] 
[TUl [2"3"|) and we obtain stochastic differential equations as limits when the time step r goes to 
0. 

Let us present the asymptotic assumption for the interaction with r = 1/n. In terms of 
the parameter n we can write the unitary operator U as 

tti \ ( 1 u \ ( L oo(n) L w (n) \ 

U(n) = exp -i-H tot = ) ' } [ } . (24) 

V n J V L on n ) L n(n) J 

Let us recall that the discrete dynamic of quantum repeated interactions is given by T4 = 
Uk ■ ■ ■ U\. In Oil], it is shown that the asymptotic of the coefficients Lij(n) must be properly 
rescaled in order to obtain a non-trivial limit for Vr nf j. With proper rescaling, is it shown in 
these references that the operator Vj nt i converges when n goes to infinity to an operator Vt 
which satisfies a Quantum Langevin Equation. When translated in our context of a two-level 
atom in contact with a spin chain, we put 

L 00 (n) = l+~W + o(~ 
n \n 

L i(n) = -^S+o (-L 
Lio(n) = —=T + o{ — 



Ln(n) = l+-Z + o(-\. (25) 
n \ n ' 
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In terms of total Hamiltonian, it is shown in [3] that typical Hamiltonian H to t which gives 



rise to such asymptotic assumption can 


be described 


as 


H tot = H Q ® I + I® ( ^ " 


+ vM 




/ 

I i o 


Hence, tor the operators IV, d, j and ^ 


we get 






W = 




7oI- 


\c*c 


z = 


-JTo- 


71 1 - 


±cc* 


s = 




-t'C 






(26) 

In the rest of the paper, we shall write all the results in terms of the operators Hq and C . 

Now, we are in position to investigate the asymptotic behavior of the different equations 
(1131 rroi [mi I23[) and to introduce the continuous models. The mathematical arguments used 
to obtain the continuous models are developed in Section 5. Before presenting the main result 
concerning the model with measurement, we treat the simpler model obtained by considering 
the limit n goes to infinity in the equation ([5]) of Proposition [TJ 

2.1 Continuous Quantum Repeated Interactions without Measurement 

In this section, by applying the asymptotic assumption, we show that the limit evolution 
obtained from the quantum repeated interactions model is a Lindblad evolution (also called 
Markovian evolution [28]). This result has been stated and proved in [3]. We recall it here 
since the more general situations treated in the current work build upon these considerations. 
The discrete Master equation ([5|) of Proposition [1] in our context is expressed as follows 

Pk+i = L(pk) = p(LooPkLoo + Li O p k L* ) + (1 - p) (L iPfc^oi + L nPkL\i) ■ (27) 
Plugging in the asymptotic assumptions (j25H . we get (here, n is a parameter) 



1 

Pfc+i = PfcH — 



n 



P [ -i[H , P ] - \{C*C, p} + CpcA + (1 - p) f-i[H ,p\ - \{CC\ p} + C*pC ) + :■( I ) 



(28) 

where [X, Y] = XY — YX and {X, Y} = XY + YX are the usual commutator and anti- 
commutator. The following theorem is obtained by taking the limit n — > oo in the previous 
equation. 

Theorem 1 (Limit Model for Quantum Repeated Interactions without Measure- 
ment) Let (p\ n t]) be the family of states defined from the sequence (p&) describing quantum 
repeated interactions. We have 

lim \\P[nt\ ~ Pt || = 0, 
where (pt) is the solution of the Master equation 

dp t = C(pt)dt, 

with the Lindblad operator C given by 

C(p) = p^-ilH ,p}-±{C*C,p} + CpC^+(l-p)(-zlH ,p]-±{CC\p} + C* P C 

(29) 
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The operator C appearing in equation ([29]) is the usual Lindblad operator describing the 
evolution of a system in contact with a heat bath at positive temperature T [3J; let us recall 
that the parameter p can be expressed in terms of the temperature T as in equation ([6]). 



2.2 Continuous Quantum Repeated Interactions with Measurement 

In this section, we present the different continuous models obtained as limits of discrete 
quantum repeated measurement models described in the equations (fT3l [To] [T9 ] 1231) , 

Although continuous quantum trajectories have been extensively studied by the second 
author in |25[ 1261 [27] , the result concerning the combination of the two kinds of measurement 
is new and the stochastic differential equations appearing at the limit have, to our knowl- 
edge, never been considered in the literature. The comparison between the different limiting 
behaviors is particularly interesting and will be discussed in detail in Section [3] 

2.2.1 The "Random environment" setup 

In Section 11.2.1] we have seen that the evolution of the little system in presence of measure- 
ment before each interaction is described by the following equation: 



Pk+i = L(Pk) + y - Vpi 1 ~ p) { L oo Pk Lq + L io Pk L\ ) 

+Vp( 1 ~ P)[Lii Pk L* n + Lqi Pk L* 01 ))x k+ i. (30) 



Using the asymptotic condition for the operator Ly(n), we get the following expression 

p fc+ i = p fc + i(£(p & ) + °(l)) +i(/C0> fc ) + o(l))x fc+1> (31) 

where the expression of C is the same as Theorem 1. The accurate expression of K is not 
necessary because these terms disappears at the limit. From the equation (|31|) . we want 
to derive a discrete stochastic differential equation. To this end, we define the following 
stochastic processes 

\ t\ 1 [ " <M 

Pn(t) = P[nt] , Fn(t) = — , W n (t) = -= X k+1 (32) 



n \ n 

k=0 



Next, by writing 

[ni]-l 

P[nt] = Po + 5^ Pk+l ~ Pk 

k=0 

and by using the equation (f3T|) and the definition ([32]) of stochastic processes, we can write 

Pn(t) = Po + [ C{ Pn {s-))dV n {s) + [ -^=£{ Pn {s-))dW n {s)+e n {t), (33) 
Jo Jo V n 

where e n (t) regroups all the o(-) terms. The equation (|33]l appears then as a discrete stochastic 
differential equation whose solution is the process (p n (t))t- 

In order to obtain the final convergence result, we shall use the following proposition 
concerning the limit behavior of the process (W n (t)). 
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Proposition 3 Let (W n (i)) be the process defined by the formula We have the following 

convergence result 

W n (t) W t , 

where =>■ denotes the convergence in distribution for stochastic processes and (Wt)t>o is a 
standard Brownian motion. 

Proof: In this case, the random variables (X^ +1 ) are independent and identically distributed. 
Furthermore they are centered and reduced. As a consequence, the convergence result is just 
an application of the Donsker Theorem |29t [3"3"] [M] . □ 

Using Proposition El we can now take the limit n — > oo in equation (j33[) . 

Theorem 2 (Limit Model for Random Environment) The stochastic process (p n (t)) ; 
describing the evolution of the small system in contact with a random environment, converges 
in distribution to the solution of the Master equation 

dp t = C(p t )dt, 

where Lindblad generator C is given in equation A29\). 

This theorem is a straightforward application of a well-known theorem of Kurtz and 
Protter [66\ [3"5] concerning the convergence of stochastic differential equations. Without the 
term l/y/n, the process (W n (t)) converges to a Brownian motion and thus the equation ([33]) 
converges to a diffusive stochastic differential equation. As the term 1/y/n converges to zero, 
this implies the random diffusive part disappears when we consider the limit. The fact that 
we recover the deterministic Lindblad evolution for a heat bath will be discussed in Section 
3. 

The next subsection contains the description of the continuous model when a measurement 
is performed after each interaction. 

2.2.2 Usual indirect Quantum Measurement 

In |25] 126] . it is shown that discrete quantum trajectories for p = converge (when n goes 
to infinity) to solutions of classical stochastic Master equations ([H [2]). These models are at 
zero temperature. The result for positive temperature (0 < p < 1) is treated in [6j. In this 
section, we just recall the result of [6] corresponding to the limit models obtained from the 
equation ([23]) . 

As it is mentioned in Section 11.2. H the final stochastic differential equations depend on 
the form of the observable. 

1. If B = aoQo + aiQi is a diagonal observable, with Qi = {q k i)o<k,i<i, we have gg = 
(/qq = 1 and all the other coefficients are equal to zero. Hence, we obtain the following 
asymptotic expression for the equation ([23]) 

p k+ i-p k = -(c(p k+ i)+o(l)^ +I(V(p fc ) + °(l))A fc+1 . (34) 

For the random variables {X k ), we have 

P + + °(1)) 

(35) 

l-p+l(g(p k ) + o{\)) 



Xk +^ = -\Z iHII with Probability 
Xk +^ = TlSiit with probability 
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In (|34p and (|35p . the exact expressions of J\f, h and g are not necessary for the final 
result. The expression of C corresponds to the Lindblad operator of Proposition 1. 

2. The other case concerns an observable B which is not diagonal. We have then < 
(/q < 1 and < < 1. The final result is essentially the same for all non diagonal 
observables B. Hence, we just focus on the symmetric case where B is of the form 



.'1/2 l/2\ 
B = a ^l/2 l/2 )+ai 



1/2 -1/2 
-1/2 1/2 

Thus, in asymptotic form, the equation (I23p becomes 

Pk+i - Pk = - (c(pk) + °(i)) + — (g(pk) + °(i))-Xfc+ij 

where Q is defined on the set of states by 



G(p) 



(36) 



- (p(Cp + P C*) + (1 - p)(C*p + pC) 
+ Tr \p(Cp + P C*) + (1 - p)(C*p + pC) 



The random variables -X^+i evolve as 



Tr^i(p. 



Tr[^ (p. 



[ X k+1 (l) = ^^0^ with probability ± + -^(m( Pfc ) 

Again, the exact expressions of / and m are not necessary for the final result. 
We define 



*j| with probability I + 



(37) 



[ni] 



[nt]-l 



Pn(t) = P[nt]> ^jnt] = — > Wn\ 



fe=0 



Depending on which type of observable we consider, we obtain two different discrete stochastic 
differential equations. 

1. In the case of a diagonal observable we have 

[nt]-l 



Pn{t) = p + Y, + o(l)) + ^(AT(ft) + o(l))^I w ) 

^ vn v n v n 7 

= p + [ C( Pn (s-)dV n (s) + / -^Af(p n (s-))dW n (s) + e n (£). 
Jo Jo V" 



2. In the same way, in the non diagonal case we obtain 

[nt]-l 

Pn(t) = PO 



l — n V 



A;=0 
t 



p + / C( Pn (s-)dV n (s) + / Q{p n {s-))dW n {s) + e n {t). 
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In these equations the terms e n (t) regroup the o(-) terms. The final results are gathered in 
the following theorem (see [2] for a complete proof). 

Theorem 3 ( Limit Model for Usual Indirect Quantum Measurement ) 

Let B be a diagonal observable. Let (p n (t)) be the stochastic process defined from the dis- 
crete quantum trajectory describing the quantum repeated measurement of B. This stochastic 
process converges in distribution to the solution of the master equation 



Let B be a non-diagonal observable. Let (p n (t)) be the stochastic process defined from the 
discrete quantum trajectory describing the quantum repeated measurement of B. This stochas- 
tic process converges in distribution to the solution of the stochastic differential equation 



where (Wt) is a standard Brownian motion. 

It is important to notice that for a diagonal observable we end up with a Master equation 
without random terms. In [26], at zero temperature, it is shown that the limit evolution is 
described by a jump stochastic differential equation. Similar evolutions for diagonal observ- 
ables will be recovered when we consider both measurements. The discussion in Section [3] 
will turn around such results. 

2.2.3 Continuous Model of Usual Indirect Quantum Measurement in Random 
Environment 

This section contains the main result of the article. To our knowledge, a random environment 
model has never been considered before in the setup of indirect quantum measurement (neither 
in discrete, nor in the continuous case). 

We treat separately the case of a diagonal observable and a non diagonal observable. We 
show that for a diagonal observable, we recover a evolution including jump random times. 
The limit evolution is although different as the case of |26j . 

Let us start with the non diagonal case. As in Section [2.2.21 we focus on the case 



dp t = C(p t )dt. 



dp t = C(p t )dt + G{pt)dW t 




In this situation, the asymptotic form of the equation (I16p 



is given by 




fc+i 



(38) 
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From the equation (|38|) . we want to derive a discrete stochastic differential equation. To this 
aim, we define the processes 

Pn(t) = P[ n t], 

r .1 1 [«*]-! 

= lj , w n (t) =-= x: 

r» * /r) ^ — * 



v fc=0 

[nt]-l [nt]-l 



and the operators 

Q(p) = v/1 - (1 - pf (Cp + pC* - Tr[p(C + C*)]p) (39) 



W(p) = s/iT^(c*p + pC-Tr[p(C* + C))p). (40) 
This way, the process (p n (t)) satisfies the following discrete stochastic differential equation 

Pn(t) = f C(p n (s-)dV n (s) + f ^=K{ Pn {s-))dW n {s) 
Jo Jo y/n 

ft 

Q{p n {s-))dW°(s) + / W{p n {s-))dWl(s) + e n {t). 
Jo 

Heuristically, if we assume that 

(W n (t),W°(t),W^t)) (W t ,W t \W?), 

where the processes iWt) and (Wj ) and (Wj?) are independent Brownian motions, the follow- 
ing theorem becomes natural (the rigorous proof is presented in Section 0J) . 

Theorem 4 (Limit Model for Indirect Quantum Measurement of non-diagonal 
observables in Random environment) Let (p n (t)) be the stochastic process defined from 
the discrete quantum trajectory (pk) which describes the repeated measurement of a non- 
diagonal observable in random environment. Then the process (p n (t)) converges in distribution 
to the solution of the stochastic differential equation 



p t = po + / C(p s )ds + [ Q{ Ps )dWl + [ W{p s )dWl (41) 
Jo Jo Jo 

where (W^) and iWf) are two independent Brownian motions. 

It is important to notice that we get two Brownian motion at the limit whereas in Theorem 
Othere is only one Brownian motion. We have already described a situation where the random 
noise disappears. 



Let us now deal with the diagonal case. In asymptotic form, the equation (fT3|) becomes 
Pk+i = Pk + -(c(pk) + + -(£(pk) + X k+X 

+ L^ fcC li -p* + f !oi - 1 ( Tr [erne*] + o(i; 



Tr [Cp k C* 

c* Pk c 

Tr [C*p k C] 



n 



+ I rSZfrt - P* + °(1) ) (lio - ^ ( H [C*PkC] + o(l))) . (42) 
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Such an equation can be written in the following way 

1 

n 



Pk+i = p k + -[C(p k )+p(-Cp k C* + Tr[Cp k C*]p k ) 

+(1 - p) ( - C*p k C + Tr [C*p k C] Pk ) + o(l 



+-(£(p k )+o(l))X k+1 



+ 1 



c Pk c* 



Tr [Cp k C* 



Pk + o(l) 1 01 + 



c* Pk c 

Tr [C*p k C] 



Pfe + o(l) 1 
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In order to define the discrete stochastic differential equation, we need to introduce the 
operator 

T(p) =C{p)+p{- CpC* + Tr [CpC*] p ) + (l-p)[- C*pC + Tr [C*pC] p) 
and the following processes 

P[nt] 



Pn{t) 
V n (t) 



^ fro 

[nt]-l [nt]~l 

£ l01 ^(i) = E ^ 
fc=0 fc=0 



(43) 
(44) 

(45) 



We obtain a discrete stochastic differential equation 

Pn(t) = Po + [ T( Pn (s-)dV n (s) + / -^£(p ri (s-))dW ri (s) 
Jo Jo V™ 

I c* Pk c 



c Pk c* 

Tr[Cp k C* 



Pk )N^s) + 



o \Tr[C* Pk C] 



p k \N 2 n (s). (46) 



Let us motivate briefly what follows concerning the convergence of (iV~(i)) and (iV^(t)) (this 
will be rigorously justified in Section 3). Let us deal with (N^(t)) for example. By definition 
of Iqi, we have 



l i(u k+1 ,ip k+1 ) = 1 with probability ± (pTi[Cp k C*} + o(l)) 
l i(wfc+i, <Pk+i) = with probability 1 - ^(pTr[Cp k C*} + o(l)) 



(47) 



Hence, for a large n, the random variable loi takes the value 1 with a low probability and 
with a high probability. This behavior is typical of the classical Poisson process [38^ [37] . 
Heuristically we can consider a counting process (iV t ) as the continuous limit of (iV*(i)). 
Since a counting process is entirely determined by its intensity ([371 EI])) we can guess its 
intensity by computing E[jV*(f)]. We have 



[nt]-l 



fc=0 

t 

E 



pTr[Cp fc C*] + o(l) 



pTr[Cp n ( S -)<T 



<%(s) + e n (*) 



(48) 
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Assuming that the processes (N^(t)) and (p n (t)) converge, we get 



E[N}]= Te 
J o 



pTi [Cps-C* 



ds. 



We thus define the limit process {N t ) as a counting process with stochastic intensity t — > 
f pTr[Cp s -C*]ds. In the same way, we assume that (N 2 (t)) converges to a counting process 
(N 2 ) with stochastic intensity t — > J" *(l — p) Tr[C* p s -C]ds. 
The limit stochastic differential equation would then be 

From a mathematical point of view, the way of defining this equation is not absolutely rigorous 
because the definition of the driving processes depends on the solution (p t ) (usually, in order 
to define solutions of a stochastic differential equation, one needs to consider previously the 
driving processes). 

A rigorous way to introduce this equation consists in defining it in terms of two Poisson 
Point processes A 1 and A 2 on M 2 which are mutually independent (see [26, 32J). More 
precisely, we consider the stochastic differential equation 

Pt = Po + J Q T {ps~)ds + J q J r i^^c~cF\ ~ Ps ~) ^^p^lCps-C^^ids^x) 

+ LL {nc%] c c} - /,s -) io<i<MTr[cift - c,Jv2(ds,&) (50) 

This allows to write the equation in an intrinsic way and, if (150)) admits a solution, we can 
define the processes 

Nt = / / ^o<x<pTr[Cp s -C*] Nl (ds,dx) and 
Jo Jm. 

= / / lo <x <(i-p)Tr[c*p^C]N 2 (ds,dx). (51) 



We can now state the convergence theorem in this context. 

Theorem 5 (Limit Model for Indirect Quantum Measurement of diagonal observ- 
ables in Random environment) 

Let A 1 and N 2 be two independent Poisson point processes on M 2 defined on a probability 
space (Q, J 7 , P). Let (p n (t)) be the process defined from the discrete quantum trajectory {pk) 
which describes the measurement of a diagonal observable A in a random environment. The 
stochastic process (p n (t)) converges in distribution to the solution of the stochastic differential 
equation 

Pt = Po + J q T iPs-)ds + jf jf \ 7^^—^ ~ Ps-^j 1 o< a; <pTr[Cp s _c*]A 1 ((is,dx) 

+ So i ( Tr C [C*p s C C] ~ Ps -) W(1 ' p) ^c^ C] N\ds, dx) (52) 
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Table 1: Different models and the corresponding continuous behavior 





No measurement 


Before 


After 


Before & After 


B diagonal 


T = Ci 
T > C p 


T = Ci 
T > C p 


T = IJ 

T > £ p 


T = IJ 
T > 2J 


B non-diagonal 


r = o id 

T > ID 


T = ID 
T > 2D 



Remark 5 It is not obvious that the stochastic differential equations Q [5|) admit a unique 
solution ( even in the diffusive case, the coefficients are not Lipschitz and the jump term can 
vanish). The uniqueness questions is treated in f25[ [2E [ [3^ 

Let us stress at this point that in this article we have focused on the particular case 
£ = C 2 . This case allows to consider observables with two different eigenvalues. In [271 E], 
situations with more than two eigenvalues are considered but only when measurements are 
performed after the interactions. The statistical model (Random environment) is not treated. 
In this article, our aim was to compare the situation with and without measurement before the 
interaction in order to emphasize the situations appearing in the case of random environment. 
The situation that we have treated is sufficiently insightful to point out the differences between 
the statistical model and the Gibbs model. Higher dimension can easily be treated by adapting 
the presentation of this article and the results of \27\ [6]; the continuous evolutions involve 
mixing between jump and diffusion evolution (see also [HI 1221 for other references on such 
types of equations). 

In the following section, we compare the different continuous stochastic Master equations 
in the different model of environment. 



3 Discussion 

The different models we have considered and the limiting continuous equations that govern 
the dynamics are summed up in Table HJ Each cell of the table contains the type of evolution 
equation in the zero temperature case (T = 0) and in the positive temperature case (T > 0). 
Hence, in what follows, the parameter p, until now supposed constant, will be allowed to vary. 
Continuous, Master equations evolutions are denoted by C p where p is the parameter related 
to the temperature (T = corresponds to p = 1). In these terms the two differential equation 
at T = are given by 

dp t =£ 1 ( Pt -)dt+ -Pt-^J (dN t -Tr[C Pt -C*)dt), (53) 

where (Nt) is a counting process with stochastic intensity J * Tr [Cp s -C*~\ds and 

dp t = Cx(p t -)dt + (Cpt + p t C* - Tr \ Pt (C + C*)l Pt )dW t , (54) 



where (Wt) is a Brownian motion. 

Note that when no measurement is performed after the interaction (the "No measure- 
ment" and "Before" columns), the type of the observable B is irrelevant. Moreover, at zero 
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temperature, the measurement before the interaction is irrelevant, since the state of the sys- 
tem to be measured is an eigenstate of the observable. Hence the last two columns contain 
identical information in the case T = 0. 

The discussion that follows is meant to provide insight about this table and on the dif- 
ferent limit behaviors that appear. We shall try, as much as possible, to provide physical 
explanations for the similarities and differences between the different models treated in the 
present work. 

3.1 Gibbs vs. Statistical models at near zero temperatures 

In order to emphasize the differences between the statistical model and the Gibbs model, we 
investigate the stochastic equations when the parameter p goes to 1, that is the temperature 
goes to zero (this fact is related to the assumption 70 > 71 in the description of the free 
Hamiltonian of £ ). In particular, we show that we can recover the zero temperature case 
from the statistical model by considering the limit p goes to 1, while it is not the case in the 
Gibbs model. This can be seen in the case of a diagonal observable. At zero temperature for 
a diagonal observable, the continuous model is given by the jump equation (j53l) . In the Gibbs 
model, for a diagonal observable, we get only the master equation dpt = C p {pt)dt. It is then 
obvious that we do not recover the equation (|53[) when we consider the limit p goes to one. 
Concerning the statistical model, i.e random environment, the limit equation is given by 



p = 1 (this result can be rigorously proved by considering the limit p = 1 in the Markov 
generator see Section HJ . Let us notice that the limit p = 1 in the diffusive evolution allows to 
recover the model at zero temperature for the diffusive evolution in both models (statistical 
and Gibbs). 

3.2 Gibbs vs. Statistical Models: absorption and emission interpretation 

In the preceding section, we have seen that the Gibbs model and the Statistical model give 
rather different continuous evolution equations, especially in the case where a diagonal ob- 
servable is measured. We are now going to provide a more complete interpretation of the 
Table [TJ To this end, we shall concentrate on the special case where TC = £ and 



This particular choice for the Hamiltonians is known as the dipole type interaction model and 
it has the property that the interaction between the small system and each copy of the chain 
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is symmetric. This will allow us to give an interpretation of the evolution of the small system 
in terms of emisions and absorption of photons. In such a setup, we shall clearly identify and 
explain the differences between the two models (Gibbs and Statistical). 

Let us start by commenting on the similarities between these models. If no measure is 
performed after each interaction, we have seen that the limit evolution is the same in both 
models. In particular, the randomness generated by the measure in the Statistical model 
disappears at the limit and we get a classical Master equation. 

The models become different when one considers a measurement, after each interaction. 
As in the previous section, the differences are more significant in the case where the measured 
observable B is diagonal. In order to illustrate the differences between the two models, we 
start by describing the trajectory of the solutions of the jump equations and by explaining 
the apparition of jumps. 

At zero temperature, the evolution equation (|53p can be re-written as 

d Pt = SxOh-W + ^[Cp^C*] ~ Pt-) dNt (56) 

by regrouping the dt terms. The solution of such a stochastic differential equation can be 
described in following manner. Let (T n ) n the jump times of the counting process (Nt), that 
is T n = 'm.i{t/N t = n}. We have then 

'< - i + 1 (4^r - ir ' s " <57) 

This expression is rigorously justified in [26] . What this means is that in the time intervals 
between the jumps, the solution satisfies the ordinary differential equation dpt = S\(pt-)dt 
and at jump times its discontinuity is given by 

Cp n -C* C Pn -C* 
PTk - PTk ~ + Tr [Cp Tk „C*] ~ PTk - ~ Tr [C m .C*] ' m 



In a similar fashion, the solution of equation (|55p satisfies 

rt ~ / Cp T i_C* \ 

Pt = J o S P (Ps-)ds + g ^ - p n -j l n <t 

c* Pn _c \ 



t- L n<t 



where for i = 0,1, the terms (Tfc) correspond to the jump times of the processes (iV t *). 
Depending on the type of the jump, the discontinuity of the solution is given by 

C prpi_C* C*prp2_C 

P n = Tr [Cp Tl _C*] ° r ^ = Tr [C*^_C] ' m 



Remark 6 Since the two Poisson point processes N 1 and N 2 are independent, on the prob- 
ability space (Q, J-, P) supporting these two processes, we have 

P [{u G n | 3k e N,T fe V) = Tl{u)}} = 0. 
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Figure 1: Experimental setup 

This means that a jump of type 1 cannot occur at the same time as a jump of type 2. We 
shall see later on that this condition is also physically relevant. 

An explicit computation with the particular value of C we considered gives 

c P c* /l 0\ C*pC ( \ 

for all states p. 

In the setup with two possible jumps, depending on the type of jump, the state of the 
small system after the jump is either ground state or the excited state. This has a clear 
interpretation in terms of the emission and absorption of photons. At zero temperature, it 
is well known that equation (|56|) describes an counting photon experiment \1'6\ [7] and that 
the a jump corresponds to the emission of a photon, which will be detected by the measuring 
apparatus. In the case where two jumps can occur, the same interpretation remains valid for 
type 1 jumps (emission of a photon). After such an emission, the state of the small system 
is projected on the ground state |eo)(eo|. Type 2 jumps are characterized by the fact that 
the state of the small systems jumps to the excited state |ei)(ei|; this corresponds to the 
absorption of a photon by the small system, which justifies its excitation. Note that the 
impossibility of simultaneous jumps of the two types (see the above remark) is physically 
justified by the fact that the small system can not absorb and emit a photon in the same 
time. 

This interpretation has a clear meaning in the discrete model. Let us consider the ex- 
perimental setup in Figure [TJ This setup, with two measuring apparatus, corresponds to the 
Statistical model. 

• At zero temperature, each copy of £ is in the ground state |eo)(eo|. In this case, the 
first measurement device will never click and only the result of the second apparatus 
is relevant. If, at the step k + 1, the second apparatus does not click, the state of the 
small system is given by pt+i = ^ooPfc-^oo/ ^[-^ooPfc-^ool- ^ n ^ ne asymptotic regime, we 
get pk+i = Pk + ^/nS\(pk) + °(1), which is an approximation of a continuous evolution. 
On the other hand, if the second apparatus clicks, then the evolution is given by Pk+i = 
L\QpkL\ Q / Tr[LiopkLi ] = CpkC* / Tr[C ' pkC*] + °(1), which corresponds to the emission 
of a photon. This corresponds to a jump, as indicated by the result of the measurement. 

• At positive temperature, both devices can click. If the first apparatus does not click, 
the state of £ before the interaction is |eo)(eo| and we have the same interpretation of 
the second measurement as before. In the other case, a click for the first measurement 
implies that the state of £ is |ei)(ei|. Now, the interpretation of the second measurement 
is the following. If we have a click, then the evolution is continuous, and the absence 




o o 
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Table 2: Physical interpretation of measurements 



App. A \ App. B 


No click 


Click 


No click 


Continuous 


Emission 


Click 


Absorption 


Continuous 



of a click corresponds to a jump of the form C* pkC / Tr[C* pkC] + o(l) (absorption of 
a photon). Let us stress that this corresponds to the inverse of the situation where no 
click occurs at the first measurement. The different cases are summarized in the Table 

m 

We are now in the position to explain the difference between the Gibbs and the Statistical 
models. In the Statistical model, the first measurement allows us to clearly identify if the 
small system absorbs or emits a photon. If we consider the same experiment without the first 
measurement device, we obtain the Gibbs model. In this setup, the information provided 
by the second apparatus is not sufficient to distinguish between a continuous evolution, an 
absorption or an emission. Indeed, as it has been pointed out in the above description, in 
order to have the exact variations of the state of the small system, it is necessary to know if 
the state of £ is |eo)(eo| or |ei)(ei| before the interaction. 

3.3 Unraveling 

In order to conclude the Section we shall to investigate an important physical feature 
called unraveling. This concept is related with the possibility to describe the stochastic 
master equations in terms of pure states. More precisely, an important category of stochastic 
master equations preserve the property of being valued in the set of pure states, that is if 
the initial state is pure, then, at all times, the state of the small system will continue to be 
pure. This property is of great importance for numerical simulations; indeed, less parameters 
are needed to describe a pure state than an arbitrarily density matrix (for a fT-dimensional 
Hilbert space, a pure state is "equivalent" to a vector that is we need 2K — 1 real parameters, 
whereas for a density matrix we need K 2 such real coordinates) . Since the expectation of the 
solution of a stochastic master equation reproduces the solution of the master equation, by 
taking the average of a large number of simulations of the stochastic master equations we get 
a simulation of the master equation. An important gain of simulation is obtained by the pure 
state property. This technique is called Monte Carlo Wave Function Method. 

When a stochastic Master equation preserve the property of being a pure state, it is 
said that the stochastic master equations gives an unraveling of the master equation (or 
unravels the master equation). In this setup, one can express a stochastic differential equation 
for vectors in the underlying Hilbert space. This equation is called stochastic Schrodinger 
equation. In this subsection, we want to show that the continuous models obtained from the 
limit of the repeated measurements before and after the interaction give rise to unraveling 
of the master equation for a heat bath whereas the unraveling property is not satisfied if we 
consider the measurement only after the interaction. Let us stress that at zero temperature, 
this property has already been established in |25} [26] (the author do not refer to unraveling 
but he shows that the stochastic master equations (|o^|) and ((Slj) preserve the property of 
being valued in the pure states set). 

In order to obtain the expression of the stochastic Schrodinger equation for the heat bath, 
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we show that the quantum trajectories can be expressed in terms of pure states. To this 
end, we show that for all k, there exists a norm 1 vector ip k G Ti. such that p k = IV'fcKV'fcl- 

(n) 

Next, by considering the process ip\ and the convergence when n goes to infinity, we get a 
stochastic differential equation for norm one vectors in 7i. Following the form of observables, 
we obtain two types of equations, which are equivalent of (f53|) and ([531) (the equivalence is 
characterized by the fact that a solution (ipt) of an equation for vectors allow to consider the 
process {\tpt)(fpt\) which satisfies the corresponding stochastic master equation). 

We proceed by recursion. Let suppose that there exists ip^ such that p k = |V , fc}('0fc|- Let 
Qi be one of the eigenprojectors of the observable B which is measured after the interaction. 
Since Qi is a one dimensional projector, there exists a norm 1 vector such that Qi = 
For j £ {0, 1}, the transitions between p k +i and p k are given by the non normalized operators 
= Tr £ i 1 ® Qi u (pk ® \ej)(ej\)U* I®Qi], for G {0, l} 2 and we have 

Zk+iUi) = Tr f U(\fa)(1> k \ ® \ej)(ej\)U* I® |i?t)(i?i|] 



Tr^ 

Tr^ 

Tr^ 
Tr £ 



^ L p i\ip k )(Tpk\Ll v ® |i?t)(i?i||e p )(ej||ej)(ej||e„)(e, 1 ||i?i)(i?i 



p,l,u,v 



® IWi 



(61) 



We can then define 



*£ p (e p ;$i)L pj ip k 



(62) 



\\Y, p { e P^i) L pj^k\\ m 

which describe the evolution of the wave function of TC. This equation is equivalent to 
the discrete stochastic Master equations in the sense that almost surely (with respect to P) 
IV'fcKV'fcl = Pk, for all k. Let us stress that, here, the normalizing factor || ^2 p {e p ; $i)L p j ip k \\ 
appearing in the quotient is not the probability of outcome. Indeed, the probability of out- 
come is || ^2 p (e p ; $i)L p j ipk\\ 2 ■ Now, we can investigate the continuous limit of this equation 
by applying the asymptotic assumptions described in Section [2j Depending on the form of 
the observable B, we obtain two different kinds of equations: 
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A jump equation (B diagonal) 

ft rl 



A = V>o+ / F P {i>s-)ds+ ( ( ( C ^L_ -ip s - j lo<x< Pt i s - Nl (dx,ds) 
Jo Jo Jm. Vv^s- / 

+ f [ - rps-) l 0<x< i- p) „ a _N 2 (dx,ds), (63) 

Jo JR V y/ v s- J 



where 



Fpilps-) = P 



(-Mo - ^{C*C + Ms _J) + ^ZC^j ^ 



+ (1 - p) ( -iH - ^{CC* + + 0^TC* j ((H) 



and [i s - = ( r ip s ^,C*C'ip s -) and u s - = (^_,CC*i/j s _). 

A diffusive equation (B non diagonal) 

h 



i>t = V>o+ / e P (^)ds+ / ^l-(l-p 2 ){C-K s I)i> s dW l s 
Jo Jo 

+ 1 VO^iC* - CsWsdWl (65) 
J o 



where 



e P (^) = p(-i#o-~(c*c-2«,c+A£r))v'< 

+ (1 - P) [~iHo - \{CC* - 2( S C* + C 2 /)) Vs, 
with k s = Re((^ s ,C'0s)) and Cs = Re(('0 s , C*ip s )). 



(66) 



By applying the Ito rules in stochastic calculus, we can make the following observation which 
establishes the connection between the equations for vectors and the equations for states. Let 
(tpt) be the solution of equation (f6"3"|) (respectively ((63]) ). then almost surely \ipt){ipt\ = Pt, for 
all t > 0, where (pt) is the solution of (j52[) in Theorem [5] (respectively (|41|) in Theorem 0]). 
Such considerations are the continuous equivalent of the remark following equation (j62|) . 

In other words, the description of the evolution of the system H in the setup with both 
measurements can be described in terms of pure states. A key property for unraveling is that 

dE[|^><^|] = £ p {E[Wt)m])dt, (67) 
for any solutions of (]6"3"]) or (|65p. 

4 Proofs of Theorems [4] and [5] 

The last section of the paper is devoted to showing that discrete quantum trajectories in 
random environment converge to solutions of stochastic differential equations (|41l [5"2"]) . We 
proceed in the following way. 



27 



In a first step, we justify rigorously the form of the stochastic differential equations pro- 
vided in Theorems 0] and Starting with the description of discrete quantum trajectories 
in terms of Markov chains, we can define the so called discrete Markov generators of these 
Markov chains. These generators depend naturally on the parameter n of the length of in- 
teraction. When n goes to infinity, the limit of the discrete Markov generators gives rise 
to infinitesimal generators. Next, these limit generators can be naturally associated with 
problems of martingale |31l 130] . The solution of martingale problems associated with these 
generators (see Definition Q] below) can be then expressed in terms of solutions of particular 
stochastic differential equations. We show that the appropriate equations are the same as 
the ones in Theorems H] and This justifies the heuristic presentation of (|41l 152 j) in Section 
l2~231 

This first step provides actually the convergence of finite dimensional laws of the discrete 
quantum trajectories to the continuous one. Finally, in a second step, we prove the total con- 
vergence in distribution by showing that the discrete quantum trajectories own the property 
of tightness (see (2"9"| [55]). 



4.1 Convergence of Markov Generators and Martingale problems 



Let us start by defining the infinitesimal generator of a discrete quantum trajectory. Let 
B = aoQo + a\Q\ be an observable where Qi = (ffjy). Let (/?&) be any quantum trajectory 
describing the measurement of the observable A in a random environment with initial state 
po- Using the Markov property (Proposition [2] in Section fl . 2 . 2 j) of (p&) on (0,C,P), we can 
consider the process (p n (t)) which satisfies 



Pn(t) = Pk 



P[Pn(0) = p\ 

k + r 



k 

- < t < 

n 



n 



' [pk+i e B\ p k 



l 

1 for all k 

Tl(p,B) for all Borel sets B, 



(68) 



where II is the transition function of the Markov chain (pk)- More precisely, the transition 
function IT is defined, for all Borel sets B, by 

n(p,B) = pTr[a o(p)]^ (p)(S)+pTV[g 01 (p)]^ l(p) (B) 

+(l-p)Tr[g 10 (p)}5 go{p) (B) + (1 -p)Tr[g u (p)]5 gi{p) (B), 

where, for i = 0, 1, we recall that 



Goi(p) 
hi(p) 



1oo L oopLo + q\ L 00 pL{ + q % m LiQpL* m + <f n Li pLi 
ql L 01 pL* 01 + q\ L 01 pL* u + q^LnpL^ + q^LupL^ 
Quip) , ,x_ Gii(p) 



Tr[Qoi(p)] 



and gi(p) 



Tr[g u (p)Y 



(69) 



It is worth noticing that the transition function II is defined on the set of states. The discrete 
Markov generator of the Markov process (p n (t)) is defined as 

A„.f(p) =nj (/(/x) - f(p))u(p,dp), 
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where E denotes the set of states and / is any function of class C 2 with compact support. 
The set of such functions is denoted by C 2 (E). In our situation, for all / E C 2 (E), we have 



Anf(p) = nUTr[g o(p)](/(/k)(p))-/(p))+pTr[goi(p)](/^i(p))-/(p) 



(1 - p) Tr[g w (p)] f(g (p)) - f(p) + (1 - p) Tr[g n (p)] f{ 9l {p)) - f(p) 



(70) 



Now, we can implement the asymptotic assumptions (J25J) introduced at the beginning of 
Section [2] and we can consider the limit of A n when n goes to infinity. In a similar way as 
Section T2.2.31 the result is divided into two parts depending on the form of the observable B. 

Proposition 4 Let A n be the infinitesimal generator of the discrete quantum trajectory de- 
scribing the measurement of a diagonal observable. We have for all f £ C 2 (E) 



lim sup\\A n f(p)-A : >f(p)\\ =0, 



(71) 



where A 1 is an infinitesimal generator defined, for all f G C 2 {E), by 



A j f( P ) 



D P f(C(p)) + 
C*pC 



f 



CpC* 



+ 



f 



Tr [C*pC] 



Tr[CpC 

- Hp) - D p f 



f(p) - Dpf 



CpC* 



pJ \Tr[CpC 
C*pC 



Tr [CpC* 



Tr [C*pC] 



Tr[C*pC]. 



(72) 



Let A n be the infinitesimal generator of the discrete quantum trajectory describing the 
measurement of the non-diagonal observable 



,'1/2 l/2\ / 1/2 -1/2 

B = ao( l/2 l/2J +ai Ul/2 1/2 



We have for all f G C 2 (E) 



lim sup\\A n f(p)-A d f(p)\\=0, 



(73) 



where A d is an infinitesimal generator defined, for all f 6 C 2 {E), by 

A d f(p) = D p f{C(p)) + \d 2 J(Q( p ), Q(p)) + ±D 2 f(W(p),W(p)), 
where Q and W are defined by the expressions OP)) and §JQ ■ 

We do not provide the proof of this proposition (similar computations are presented in 
great detail in [27]). Now, we can introduce the martingale problem associated with the limit 
generators of the above Propositions) To this aim, we denote {!F^ ), the filtration generated 
by a process (pt), where T% = cr(p s , s < t) for all t > 0. 
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Definition 1 Let (Q,J-,¥) be a probability space. Let i G {j, d} and let po be a state on E. 
A solution associated with the problem of martingale (A l ,po) is a process (pi) such that, for 
all f G C 2 , the process (M%(f)) defined by 

Ml(f) = f(pl)-f( Po )- f A l f(pl)ds 

Jo 

is a [Tl ) martingale with respect to P. 

Usually, solutions of stochastic differential equations are used to solve the problems of mar- 
tingale [3TJ [HUJ. In our context, we recover the stochastic differential equations (|4*Tl |5"2"|) 
introduced in Theorems H] and [5l Let us start by the non diagonal case. 

Theorem 6 (Solution of the Problem of Martingale for a Non-Diagonal Observ- 
able) Let (Q,J-,'P) be a probability space which supports two independent Brownian motions 
(Wf) and (Wt). Let A d the infinitesimal generators corresponding to the discrete quantum 
trajectory describing the measurement of 



,1/2 l/2\ / 1/2 -1/2 
A - Ao( l/2 1/2 J + Al 1-1/2 1/2 



Let (po) be any state. The solution of the problem of martingale associated to (A d ,po) is given 
by the solution of the following stochastic differential equation 

Pt = P o+ f C(Ps)ds + f Q{ Ps )dWl + f W{ Ps )dWl (74) 
Jo Jo Jo 

The equivalent theorem in the diagonal observable case is expressed as follows. 

Theorem 7 (Solution of the Problem of Martingale for a Diagonal Observable) 

Let (0,,J-,P) be a probability space which supports two independent Poisson Point Process N l 
and N 2 . Let A^ the infinitesimal generators corresponding to the discrete quantum trajectory 
describing the measurement of a diagonal observable. Let po be any state. The solution of the 
problem of martingale associated to (A 3 , po) is given by the solution of the following stochastic 
differential equation 

Pt = Po + f Q T iPs-)ds + jf ^ Tr ^~ C c *j ~ P'-j 1 o< a; <pTr[Cp s _c*]^ rl (rfs,dx) 

T^[C*p C C] ~ Ps ~) 1 °< :E <( 1 -P) Tr [ c '>^ c ] iV2 ^ s ' dx ^ ( 75 ) 

These theorems can be proved by using Ito stochastic calculus (see [27] for explicit computa- 
tions). 

In order to complete the study of the limit infinitesimal generators, we express a uniqueness 
theorem of solutions for the problems of martingale. Moreover this result is essential to prove 
the final convergence theorem. 
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Proposition 5 Let (po) be a state and let A 1 , i = 0, 1 be a generator defined in Proposition 
5. The problem of martingale (A l ,po) admits a unique solution in distribution. It means that 
two solutions of the martingale problem (A l ,po) have the same law. 

This proposition is actually a consequence of the uniqueness of solution for the stochastic 
differential equation associated with A 1 . Complete reference about Markov generators and 
problems of martingales (uniqueness, existence) can be found in [30J. 

The next section contains the final convergence result. 



4.2 Tightness Property and Convergence Result 



We prove that discrete quantum trajectories (p n (t)) have the tightness property (also called 
relative compactness for stochastic processes). Next, we show that the convergence result of 
Markov generators (Proposition U]) implies the convergence of finite dimensional laws. The 
tightness property and the finite dimensional laws convergence imply then the convergence in 
distribution for stochastic processes [29J. 

Concerning the tightness property, we have the following result. 

Proposition 6 (Tightness) Let (p n (t)) be any quantum trajectory describing the repeated 
quantum measurement of an observable A (diagonal or not). There exists some constant Z 
such that for all t± < t < t 2 



E 



\Pn(t 2 ) ~ Pn(t)\\ 2 \\Pn(t) ~ Pn(h)\\ 2 < Z(t X - t 2 ) 



As a consequence, the sequence of discrete processes (p n (t)) is tight. 



(76) 



In order to see that the property (176]) implies the tightness property, the reader can consult 
. Before to prove the Proposition El we need the following Lemma. 

Lemma 1 Let (p^) be the Markov chain describing the discrete quantum trajectory defined 
by the repeated quantum measurement of an observable A. Let 

M^ = a{ Pj ,j<r}, 

and let (r, I) E N 2 such that r < I. Then there exists a constant Ka such that 



E 



\pi -prf/M 



(n) 



<K A x 



I — r 



n 



Proof: We just treat the case where B is diagonal (similar reasoning yield the non diagonal 



case). Let us start with the term defined by E \\pi — Pr\\' 2 /-M[™\ 



We have 



E 



\PI-Prf/Mt\ 



E 



E 



E 



i 



M 



(n) 
l-l 



E \Mpi-i) - Pr\\ 2 P Tr[^(pi-i)] I MfX 
+ E Il^(^-l) " Prf (1 - V) Tt[&i(Pi-i)] / Mf\ 



(n) 



(77) 
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With the asymptotic description of <pi and 9{, we have for the first term in the right side of 
expression (|77j) 



E 

= E 
+E 

< pE 



J] \\MPl~l) ~ Pr\\ 2 pT*[Q0i(Pl-l)] I Mf\ 

i 

\\pi-i + - (A)(Pi-i) + o(l))) - p r || 2 p(l - -(Tt[C«_iC7*] + o(l))) / Mt\ 
n n I 

\\h{px-x) - Prf-mCpi-tC*] + o(l))) / Mt\ 
n I 



1 



+- x E 

n 

+— x E 

n 



Ipi-l-Prf/M™ 

1 



£o(w-i) + o(i))||Mi - £ Cft[Cp,_iC*] + o(i)) / m£\ 



n 



11 



f 2 (p { _!) - ^fc^c^c*] + o(i))) / .M^ 



As the discrete quantum trajectory (p^) takes values in the set of states which is compact 
and as the function defined on the set of state p i — > f2(p)( Tr[CpC*] + o(l)) is continuous, 
there exists a constant Zi such that, almost surely 



E 



Y,\\Hpi~i)-Pr\\ 2 v^[Q^( P i-i)] /M™ 



\PI-1-Pr\\ 2 M£\ 



+ (78) 
n 



In the same way there exists a constant Z2 such that 



E 



E l|0<(ft-i) " Pr\\ 2 pTr[g u (pi-i)} I Mf\ 



< (l-p)E 



\PI-1- Prf Mf\ 



Z2 

n 



Finally, for an appropriate constant Z, we have almost surely 



E 



\PI-Pr\\ 2 /Mt\ 



< E 



As a consequence, by remarking that 



E 



\Pl ~ Pr\ 



E 



E 



\ Pl -! - Prf / M£\ 



\ PI -Pr\\ 2 Mt\ 



z 

+ -. 

n 



(79) 



(80) 



by induction, we have 



E 



\Pl ~ Pr 



I 2 



<K A - 



l-r 



n 



In the non-diagonal case, the computation and estimation are similar and the Lemma holds. 
□ 



Proposition [6] follows from this lemma. 



32 



Proof: (Proposition [6]) Thanks to Lemma 1, for all quantum trajectories (p n (t)), we 
have: 

E[\\p n (t 2 ) - p n (t)\\ 2 \\p n (t) - p n (h)\\ 2 ] 



E 



E 



\Pn([nt 2 }) - p n ([nt})\\ 2 /M { ^ t] \ \\p n ([nt}) - p n ([nt 



< K A ([nt 2 ] - [nt]) E 



n 



E 



\ P n{[nt])- Pn {[nt x ])r/Ml tl 



< 



K A ([nt 2 ] - [nt])K A ([nt] - [nti]) 



n 



n 



< Z A (t 2 -h) 2 , 
with Z A = 4(K A ) 2 and the result follows. 



□ 



Since the tightness property holds, it remains to prove that the finite dimensional laws 
converge. This result follows from the following proposition. 

Proposition 7 Let po be a state. Let (p n (t)) be a quantum trajectory describing a repeated 
quantum measurement of an observable A. Let A 1 , i = 0, 1 be the associated Markov generator, 
we have 



lim E 

n— >oo 



f(p n (t + S ))-f(p n (t) 



t+s 



^7(Pn(s)))n^n&)) 



i=l 



(81) 



for all m > 0, for all < t\ < t 2 < • • • < t m < t < t + s, for all functions (0j)i=i ...,m an d f or 
all f in C 2 . 

Proof: Let (p n (t)) be any discrete quantum trajectory and A 1 the associated generator. Let 
(J 7 ^) denote the natural filtration of the process (p n (t)), that is 

(Ti n) )=o{p n {s),s<t} = M^ t] . 

For m>0,0<ti<t 2 <...<t m <t<t + s and f,6±, . . . ,6 m £ C 2 , we have 



E 



E 



E 



f( Pn (t + s))-f( Pn (t)- 
f{ Pn {t + s))- f{p n {t) 



A'fipnWUllOiiPniti)) 



i=l 



t+s 



A 1 f (puis)) /J? 



Y[0i( P n(k)) 



i=l 



(82) 



Let us now estimate the term E 



J{ Pn {t + s)) - f{p n {t) - J t t+S A'fipnis))) /Jf\ . To this 
end, from the definition of infinitesimal generators, we can notice that the discrete process 
defined for all n by 



fc-i 

f(p n (k/n)) - /(po) - -A^fipnU/n)) 



(83) 



3=0 



is a (^ r i ^ /ra )-martingale (this is the discrete equivalent of solutions for problems of martingale 
for discrete processes). 
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Now, assuming r/n <t< (r + l)/n and l/n < t + s < (l + l)/n, we have = F^ n . The 
random states p n (t) and p n (t + s) satisfy then p n (t) = p n (r/n) and p n [t + s) = p n {l/n). The 
martingale property (|83H implies then 



E 
E 

E 



f(p n (t + s ))-f(p n (t)r? 



[f(Pn(l/n)) - f( Pn (k/ 
1-1 

E- A nf(Pn(j/n))/^ /n 



r/n 



Tl 

j=k 
t+s 



E 

+E 



A n f( Pn (s))ds F? 



n 



AifiPnit)) + 



(t + s) I A n f(p n (t + s)) / 



(84) 



As a consequence, we have 



E 



t+s 



< E 



f( Pn (t + s))-f(p n (t)- 

^ A l J{p n {s))-A i f{ Pn {s)ds 



1=1 



n 

8=1 



+E 



< Msup 
pes 



t- A n f( Pn (t)) + ( [ ^^ -(t + s)) A'JlnJl + *)) 



^/(p)-AV(p) 



n 

i=l 



L 


Afip) 


H sup 


« pes 





(85) 



where M and L are constants depending on \\hi\\ and s. Thanks to the condition of uniform 
convergence from Proposition 2, we obtain 



lim 

n— >oo 



E 



f( Pn (t + S)) - f( Pn (t) 



t+s 



Aif(p n (s))) f[ei(p n (U)) 



1=1 



0. 



(86) 



□ 

We finish by showing that Propositions [6] and [7] implies the convergence in distribution. 
Indeed the tightness property, which is equivalent to relative compactness for the Topology of 
Skorohod [33l[29], implies that all converging subsequence of (p n (t)) converges in distibution 
to the solution of the problem martingale (A l ,po). In other terms, let (Yt) be a limit process 
of a subsequence of (p n (t)), Proposition [7] implies that 



E 



f(Y t+s ) - f(Y t ) 



t + s 



A l f(Y s )ds)\{e i {Y ti 



i=l 



0. 



(87) 



for all m > 0, for all < t\ < t% < . . . < t m < t < t + s, for all functions (0i)i=i,..., m and for 
all / in C^. As a consequence (Yt) is a Markov process (with respect to its natural filtration 
(J-J)), which is also a solution of the martingale problem (A l ,po). Now, the uniqueness of 
the solution of the problem of martingale (Proposition [S]) allows to conclude that the discrete 
quantum trajectory converges in distribution to the solution of the problem of martingale. 
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